COLOUR CONFINEMENT IN THE 
LATTICE LANDAU GAUGE QCD SIMULATION 



Sadataka FURUI 
School of Science and Engineering, Teikyo University, 
E-mail: furui@liberal. umb. teikyo-u. ac.jp 

Hideo NAKAJIMA 
Department of Information Science, Utsunomiya University, 
E-mail: nakajima@is. utsunomiya-u. ac.jp 

The colour confinement criterion proposed by Kugo and Ojima is tested in the 

lattice Landau gauge QCD simulation. The renormalization effects are studied 

by measuring the gluon propagator, ghost propagator, three gluon vertex and the 

ghost-antighost-gluon vertex. The running coupling a s from ghost-antighost-gluon 

vertex in the infrared region yields A ~ IGeV, consistent to that from three 

° J MOM 
gluon vertex in high momentum region. 



1 The confinement signal and the definition of the gauge field 

Two decades ago, Gribov pointed out a possible mechanism of colour confine- 
ment in Coulomb gauge or Landau gauge QCD via infrared divergence of the 
Faddeev-Popov ghost propagatom. At nearly the same time, Kugo and Ojima 
proposed a criterion for the absence of coloured massless asymptoptic states 
in Landau gauge QCD using the BRST symmetrp. 

We study the confinement signal in the gluon propagator, the ghost prop- 
agator and the Kugo-Ojima parameter, and their dependence on the gauge 
field. Usually, the gauge field on lattice A^x) is defined from the link variable 
Up(x) as A x ^ = \(JJ x ,n - Ul ^)\traceiess part-, which we call [/-linear version. 
A more natural definition^ is U x ^ = expA x ^, A x „ = — A x .^, which we call 
logf/ version. 

We observe that the Kugo-Ojima parameter u ab (0) at (3 = 5.5 is about 
-0.7 in the log U and about -0.6 in the [/-linear version. 

If the configuration is in the core regioru, the tensor G°^ v = tr(\ a D^ _g D {—D v )\ h ) xy 
divided by N 2 — 1, where N is the number of colours—is expected to ap- 
proach a function E(U) defined by the optimizing function! They are E{U) = 

J2i -k^ e trt/; in [/-linear version, and N ±_ 1 J2i a tr (^ at S(Ai)X a ^j , where Ai = 
adjAi and S(x) — th ^/ 2 ) » m 1°§^ version. 

In the table below ei and e2 stand for e = ( E \/' ) , in our 16 4 lattice 
simulation of the [/-linear and the log U version of the gauge fields, respectively. 
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Figure 1: The dependence of the Kugo- 
Ojima parameter measured along x,y,z and 
t axis resp. on the definition of the gauge 
field. j3 = 6.0, 16 4 . The triangles are log U, 
and diamonds are U— linear version. 



Figure 2: The ghost propagator as a func- 
tion of the lattice momentum. /3 = 6.0, 16 4 . 
Triangles are log U definition and diamonds 
are [/-linear definition. The dashed curve is 
1.162/p 2545 



We define the horizon function h = ( H y ) / (3(iV 2 — 1)) = c — e/d. In general 
we expect h < and in the continuum limit, we expect h — when the 
configuration is in the core region. 



Table 1: fi dependence of the Kugo-Ojima parameter c, trace e divided by the dimension d, 
and h = c — e/d. The suffix 1 corresponds to the u-linear and 2 corresponds to the log U 
version. Data are those of 16 4 , except (3 = 5.5 [/-linear data, which are those of 8 4 . 



a 


Cl 


ei/d 


h 


C2 


e 2 /d 


h 


5.5 


0.570(58) 


0.780(3) 


-0.21 


0.712(18) 


0.908(1) 


-0.20 


6.0 


0.576(79) 


0.860(1) 


-0.28 


0.628(94) 


0.943(1) 


-0.32 



The gluon propagator is infrared finite and the absolute value in log U 
version is about 20% larger than that in {/-linear version The corresponding 
difference in the ghost propagator is about 10%. 



2 The QCD running coupling 

The QCD running coupling a(/i) = g 2 /47r can be measured from the three 

n q (3)/ 2 2 2 )z 3 ^ 2 ( 2 ) 

gluon verticeli as, g^ 2 ) = Ga m {p ^ G ^2) 2 ) £ A 2 j ■ Since the lattice data 

of Ga^ is infrared finite, An contrast to the conjecture that the gluon propa- 
gator is infrared vanishing], g(n 2 ) decreases as fi decreases as 1.5/i. This be- 
haviour does not agree with the results of Dyson-Schwinger approaclfl, which 
suggest that a s {n) monotonically increases to a finite constant as fi goes to 0. 
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Figure 3: The gluon propagator as a func- 
tion of the lattice momentum. j3 = 6.0, 16 4 . 
Triangles are log U version and diamonds 
are [/-linear version. 50 samples. 



Figure 4: The QCD running coupling con- 
stant as a function of momentum fi. Trian- 
gles and diamonds are g 2 /4n, the same as 
Figpj, stars and boxes are g 2 /4n and g 2 /4ir, 
respectively of log U smeared version. 



The running coupling can be obtained also from the ghost-antighost-gluon 

r ~t 2\ G, (3) (p 2 , Pf 2 ,p a 2 )z}. /2 (p, 2 )Z 3 ( fJ , 2 ) ,. . . , 

coupling g(n') = ^ {2 ) {p . 2)Gg w (p/ 2 )Ga (2) (pc2) ■ The preliminary results of 

g 2 /An and g 2 /An at symmetric momentum points (p^) = ~ sin(n M 7r/i), n 2 ^ = 
10 have a peak which depend on definition of the but at other momentum 
points g 2 /Ait and g 2 /An are not so different. When the physical scale is fixed 
by a' 1 = 1.91 ± O.lGeV, A Jl ~ / calculated from the g 2 /Air at J2 M n fi = 8 and 
16 are about lGeV, which are ccuasistent to that obtained from the three gluon 
vertex in high momentum regioiH. 
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